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Abstract 

We prove a law of large numbers and a functional central limit 
theorem for multivariate Hawkes processes observed over a time in- 
terval [0,T] in the limit T — > oo. We further exhibit the asymptotic 
behaviour of the covariation of the increments of the components of 
a multivariate Hawkes process, when the observations are imposed by 
a discrete scheme with mesh A over [0, T] up to some further time 
shift r. The behaviour of this functional depends on the relative size 
of A and r with respect to T and enables to give a full account of 
the second-order structure. As an application, we develop our results 
in the context of financial statistics. We introduced in [17] a micro- 
scopic stochastic model for the variations of a multivariate financial 
asset, based on Hawkes processes and that is confined to live on a tick 
grid. We derive and characterise the exact macroscopic diffusion limit 
of this model and show in particular its ability to reproduce important 
empirical stylised fact such as the Epps effect and the lead-lag effect. 
Moreover, our approach enable to track these effects across scales in 
rigorous mathematical terms. 
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1 Introduction 

1.1 Motivation and setting 

Point processes have long served as a representative model for event-time 
based stochastic phenomena that evolve in continuous time. A compre- 
hensive mathematical development of the theory of point processes can be 
found in the celebrated textbook of Daley and Vere- Jones [8], see also the 
references therein. In this context, mutually exciting processes form a spe- 
cific but quite important class of point processes that are mathematically 
tractable and widely used in practice. They were first described by Hawkes 
in 1971 [13, 12] and according to [8]: "Hawkes processes figure widely in 
applications of point processes to seismology, neurophysiology, epidemiol- 
ogy, and reliability... One reason for their versatility and popularity is that 
they combine in one model both a cluster process representation and a sim- 
ple conditional intensity representation, which is moreover linear. It comes 
closest to fulfilling for point processes, the kind of role that the autoregres- 
sive model plays for conventional time series." Informally, to a multivariate 
d-dimensional counting process N = (N\, . . . , N^) with values in N d is asso- 
ciated an intensity function (Ai, . . . , Xd) defined by 



where P stands for probability and ^ is the sigma-field generated by N up 
to present time t. A multivariate Hawkes process has intensity 



and is specified by fii S M+ = [0, oo) and for i = l,...,d, the (pij are 
functions from R + to K + . More in Section 2 below for rigorous definitions. 

The properties of Hawkes processes are fairly well known: from a prob- 
abilistic point a view, the aforementioned book of Daley and Vere-Jones [8] 
gives a concise synthesis of earlier results published by Hawkes [13, 12, 14] 
that focus on spectral analysis following Bartlett [4] and cluster representa- 
tion, see Hawkes and Oakes [15]. From a statistical perspective, Ogata stud- 
ied in [25] the maximum likelihood estimator and showed that parametric 
experiments generated by the observation of Hawkes processes are regular 
in the sense of Fisher information (see for instance [9] for a modern formu- 
lation of statistical regularity). Ogata also studied the numerical issue of 
exact simulation of Hawkes processes [26] , and this topic has known several 
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developments since, see [22] for a review. Recently, the nonparametric esti- 
mation of the intensity functions has been investigated by Reynaud-Bouret 
and Schbath [27] and Al Dayri et al. [3] in dimension d = 1. 

However, in all these papers and the references therein, the focus is on 
the "microscopic properties" of Hawkes processes, i.e. their infinitesimal 
evolution, possibly under a stationary regime or for a large time horizon 
[0, T] in order to guarantee a large number of jumps for statistical inference 
purposes. In the present work, we are rather interested in the "macroscopic 
properties" of Hawkes processes, in the sense of obtaining a limit behaviour 
for the multivariate process (Nt v )v£[o,i] as T — )• oo, for a suitable normal- 
isation. Our interest originates in financial data modelling: in [17], we in- 
troduced a stochastic model for the price S = (Si, . . . ,S n ) of a multivariate 
asset, based on a d = 2n dimensional Hawkes process of the form (1), with 
representation 

S 1 = N 1 - N 2 , S 2 = N 3 -N A ,..., S n = N d _! - N d . (2) 

In this context, the fact that the Si take values on Z accounts for the discret- 
ness of the price formation under a limit order book. If we take fJ, 2 i-i = l^2i 
for every 1 < i < d, the process S is centred and the mutually exciting 
properties of the intensity processes Aj under (1) allow to reproduce em- 
pirically microstructure noise and the Epps effect, as demonstrated in [17]. 
Microstructure noise - a major stylised fact in high frequency financial data 
(see e.g. [10, 23, 24, 28, 29, 21]) - is characterised by the property that in 
microscopic scales 1 , an upward change of price for the component Si is more 
likely to be followed by a downward jumps of Si and vice versa. Similarly, 
a jump of a component iVjji— i or N 2 { of Si at time t will increase the con- 
ditional probability, given the history up to time t, that a jump occurs for 
the component N 2 j-i or N 2 j of Sj through the exciting effect of the kernels 
V?2i-i,2j-i, <P2i-i,2j) ¥>2i,2j-i or (f 2 i y2 i, thus creating a certain dependence 
structure and the Epps effect, another well-document stylised fact implying 
that the empirical correlation bewteen two assets vanishes on microscopic 
scales [11]. However, both microstructure noise and the Epps effect vanish 
at coarser scales where a standard diffusion regime dominates. Thus, a cer- 
tain macroscopic stability property is desirable, namely the property that 
S behaves like a continuous diffusion on coarse scales. Also, one would like 
to be able to track the aforementioned microscopic effects produced by the 
kernels (fij in the diffusion limit of S. This is the main topic of the paper 

1 when the data are sampled every few seconds or less. 
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and also its novelty, as far as the modelling angle is concerned. 

The importance of mutually exciting processes has expanded in financial 
econometrics over the last years, see [7, 5, 20, 16] among others. Prom a 
wider perspective however, a mathematical analysis of the limiting behaviour 
of Hawkes processes has a relevance beyond finance: apart from the interest 
of such a study for its own sake, one could also presumably apply asymptotic 
results in other application fields. These include seismology of course, for 
which Hawkes processes were originally introduced (see [31] for instance) and 
also the citation we used above from the monograph [8]; in a more speculative 
or prospective way, one could mention recent related areas such as traffic 
networks [6] or genomic analysis [27] where mutually exciting processes take 
a growing importance. 

1.2 Main results and organisation of the paper 

Because of this versatility in the use of Hawkes processes, we present our 
results in a fairly general setting in Sections 2, 3 and 4 with relatively mini- 
mal assumptions. It is only in the subsequent Section 5 that we consistently 
apply our results to the multivariate price model S described by (l)-(2) for 
financial data modelling. 

In Section 2, we recall the construction of Hawkes processes in a rigourous 
and general setting. We state in particular the fact that as soon as the 
kernels (fij that define the intensity process A are locally integrable, the ex- 
istence of ./V with uniqueness in law is guaranteed on a rich enough space, 
as follows from the theory of predictable projection of integer-valued ran- 
dom measures of Jacod [18]. We state in Section 3 a law of large numbers 
(Theorem 1) and a functional central limit theorem (Theorem 2). The law 
of large numbers takes the form 

sup ||r _1 iVT„-'y(Id-K)-V|| ^° as T^oo (3) 
«e[o,i] 

almost surely and in L 2 (P). We set || • || for the Euclidean norm on W 1 . The 
limit is described by jjL = (fj,\, . . . ,/id) identified with a column vector and 
the d x d matrix K = J* °° (pit) dt, where ip = ((fij) is a matrix of functions 
from to R+. The convergence (3) holds under the assumption that the 
spectral radius of K is strictly less than 1, which guarantees in particular all 
the desirable integrability properties for N. As for a central limit theorem, 
under the additional condition that t l l 2 tp(t) is integrable componentwise, a 
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consequence of Theorem 2 is the convergence of the processes 

Vrfr^Nrv - v(id - K)- V) , ve[o, i] 

to 

(id - K)- 1 ^ 1 / 2 ^, ve[o,i\. 

in law for the Skorokod topology. The limit is described by a standard el- 
dimensional Brownian motion (W 1) )„ g [ ji and the d x d diagonal matrix £ 
defined by Sjj = ((Id — K)" 1 ^)^. The proof of Theorems 1 and 2 relies 
on the intensity theory of point processes: using that Nt — f Q X s ds is a 
c?-dimensional martingale, if we set Xt = Nt — E(N) with E(-) denoting 
expectation, we then have the fundamental representation 

X t = M t + [ ip{t-s)X s ds 
Jo 

of X as a sum of a martingale and a convolution product of <p with X itself. 
Separating the two components, we can then take advantage of the powerful 
theory of limit theorems of semimartingales, as exposed in the comprehen- 
sive book of Jacod and Shiryaev [19]. 

In Section 4, we consider yet another angle of study that is useful for 
applications. Given two square integrable point processes N and N' with 
value in N d , setting X t = N t - E(N t ) and X[ = N' t - E{N' t ), one can define 
the empirical cross-correlation of N and N' at scale A over [0, T] as 

v a>t (n, JV) = -J (XiA - * (i -i)A) {xL - x{ t _ 1)A ) T , 

i=i 

where Xt and X[ are identified as column vectors and (-) T denotes trans- 
position in K d . In rigorous terms, V^tC-^j N') shall rather be called the 
empirical correlation matrix of the increments of X and X' sampled at scale 
A over [0,T]. This object is of major importance in practice: viewed as a 
function of A, it reveals the correlation structure across scales between N 
and N'. It is crucial for the understanding of the transition of (Nt v ) V £[o,i] 
from a microscopic regime of a point process to the macroscopic behaviour of 
a Brownian diffusion, up to appropriate normalisation, as T — > oo. More pre- 
cisely, given a time shift r = ry S 1, we focus in the paper on the r-shifted 
empirical correlation matrix of N at scale A, namely 2 Vaj^JV^.). We 

2 This is actually a more general object since we can formally recover Va,t(N, N') from 
Va,t(N",N" + .), where iV" = (N,N') is a 2d-dimensional point process by letting r — > 
in V a ,t{N" , ]<'+.)■ 
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prove in Theorem 3 that Va,t(^ N t+ .) is close in L 2 (P) to a deterministic 
explicit counterpart 

v A , r = I (l- l ±^) + (ld5 {ds) + ^{s)ds)X(ld5 (dt)+tP(t) T dt) 

Jm.\ 

under appropriate conditions on r and A relative to T as T — > oo. The 
approximation va, t is described by S that already appears in Theorems 1 
and 2 and tp = Yl n >i ( Pm where <~p n is the matrix of n-fold convolution prod- 
uct of ip. This apparently cumbersome formula paves the way to explicit 
computations that yield crucial information about the dependence structure 
of functionals of N, as later developed in the paper. 

Section 5 is devoted to some applications of Theorems 1, 2 and 3 to the 
study of the multivariate price model S described by (l)-(2) for financial data 
modelling. We elaborate on particular examples the insight given by va, t - 
Concerning macroscopic limits of a univariate price S = N± — N2 obtained 
for d = 2, we derive in Proposition 1 the limiting variance of the rescaled 
process T~ 1 / 2 {N\^tv — -^2,T«) in the special case where <pi^ = </?2,i = <P 
and ipi i = (f2 t 2 = 0. This oversimplification enables to focus on the phe- 
nomenon of microstructure noise and reveals the macroscopic trace of in 
the variance of the limiting Brownian motion. We discuss the theoretical 
implications of the obtained formulas in statistical finance. In the same way, 
we explicit in Proposition 2 the macroscopic correlation structure obtained 
with a bivariate price process S = (S\, S2) = (Ni — N2, N3 — N4) when only 
cross-excitations tpi^ = and (^2,4 = ^2,4 are considered, the other com- 
ponents ifij being set to for simplicity. Again, we obtain an illuminating 
formula that characterises the influence of the microscopic dynamics on the 
macroscopic correlation of the Brownian limit. This enables in particular in 
Proposition 3 to demonstrate how Hawkes process can account for the Epps 
effect [11] and also the lead-lag effect between two financial assets that has 
been given some attention in the literature recently [30, 2, 1]. 

Section 6 develops some preliminary tools for the proofs. Section 7, 8 
and 9 are devoted to the proof of Theorems 1 , 2 and 3 respectively. Some 
technical results and computations of the application Section 5 are delayed 
until an appendix. 
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2 Multivariate Hawkes processes 



Consider a measurable space (fi, J^") on which is defined a non-decreasing 
sequence of random variables (T n ) n >x taking their values in (0, oo], and such 
that T n < T n+ \ on the event {T n < 00} , for all n > 1. Let (Z n ) n >i be a 
sequence of discrete random variables taking their values in { 1, . . . , d} for 
some positive integer d. Define, for t > 

Ni, t = ^2 1 {T n <t}n{z n =i}- 

n>l 

Remark that 2Vjo = by construction. We endow fi with the filtration 
0^t)t>o where is the <r-algebra generated by the random variables iVj s , 
s < £, 1 < z < cZ. According to Jacod [18], for any progressively measurable 
non-negative processes (Ai^)t>o, . . . , (^d,t)t>o satisfying 

Xi tS ds < 00 almost-surely, 

there exists at most one probability measure P on (fi, j^oo) such that the 
compensator (or predictable projection) of the integer-valued random mea- 
sure 

N(dt,dx) = ^2l {Tn<oo} 5 {TntZn) (dt,dx) 

n>l 

on (0, 00) x {1, . . . , d} is 

d 

u(dt, dx) = ^i,t dt (8) Si(dx), 

i=l 

where 5 is the Dirac mass. In other words, for all n > 1, for all i G {1, . . . , d}, 
the process 

rtAT„ 

N i,tAT„ ~ / ds 

Jo 

is a (j^)-martingale. This implies that the law of the (i-dimensional process 
(Aq, . . . , N,i) is characterised by (Ai, . . . , A^). Moreover if Q is rich enough 
we have the existence of such a probability measure P. 

Definition 1. We say that A = (Aq, . . . ,Nd) is a multivariate Hawkes pro- 
cess when 

~ d 

Aj,t = m + / y2 Witt ~~ s )dN j>s (4) 
where /U, £ and is a function from R + to M+. 
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We have a non-explosion criterion, as a consequence for instance of equal- 
ity (10) in Lemma 2, Section 6 below. 

Lemma 1 (Non-explosion criterion). Set = lim n T n . Assume that the 
following holds: 

j ipij(s)ds < oo for all i,j,t. (5) 
Jo 

Then = oo almost surely. 



3 Law of large numbers and functional central limit 
theorem 

On a rich enough probability space (O, j^, P), we consider a Hawkes process 
TV" = (N t )t>o according to Definition 1, satisfying (5) and specified by the 
vector 

H = (hi, ...,Hd) 
and the d x <i-matrix valued function 

Note that in this setting, we do not assume a stationary regime for N . 
Consider the assumption 

For all i,j we have J °° ipij(t)dt < oo and the spectral radius fAl) 
p(K) of the matrix K = J °° <p(t) dt satisfies p(K) < 1. 

First we have a law of large numbers in the following sense: 

Theorem 1. Assume that (Al) holds. Then N t G L 2 (P) for all t > and 

we have 

sup Hr^ 1 ^ - v (Id - K)~V|| ^° as T ^ co 

we [o,i] 

almost- surely and in L 2 (P). 

Next we have an associated functional central-limit theorem. Introduce 
the functions (p n defined on K+ and with values in the set of d x d-matrices 
with entries in [0, oo] by 

Vi = Vn+iit) = I f^- s )fn{s) ds, n > 1. (6) 
J o 
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Under (Al) we have J °° ip n {t) dt = K n hence the series Yln>l converges 
in L l (dt). We set 

^=J>n (7) 

n>l 

Theorem 2. Assume that (Al) holds. We have 

E(N t )=tn+(^J ip(t - s)sdsjfi 
where ij) £ L l (dt) is given by (7). Moreover, the processes 
■ ]=(n Tv -E(Ntv)), «e [0,1] 



converge in law for the Skorokod topology to 

(id- K)- 1 *] 1 / 2 ^, we [0,1] 

as T — > oo, where (W v ) v€ [q^ is a standard d-dimensional Brownian motion 
and S is the diagonal matrix such that Yin = ((Id — K) _1 /i)j. 

Consider now the following restriction on ip: 

/•oo 

/ fit) t 1 ^ 2 dt < oo componentwise. (A2) 
Jo 

Using Theorem 1 and Assumption (A2), we may replace T~ 1 E(Nt v ) by its 
limit in Theorem 2 and obtain the following corollary. 

Corollary 1. Assume that (Al) and (A2) hold. Then the processes 

Vt^n Tv - W (id-K)-V), [0,1] 

converge in law for the Skorokod topology to 

(Id- K)^ 1 / 2 ^, ue[0,l] 

as T — > oo. 
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4 Empirical covariation across time scales 

For two square integrable counting processes N and N' with values in N rf , 
set 

X t = N t -E(N t ), X' t = N' t -E{N' t ). 

The empirical covariation across time scales of N and N' is the process 
Va,t(N, N'), T > 0, A > 0, taking values in the set of d x d matrices and 
defined as 

j L^/Aj 

Va,t(N,N') = - (^a-X (4 _ 1)A )(X 4 ' a -X ( ' 4 _ 1)A ) t 

i=l 

where (A^a - -Xyi-i)a) is identified as a column vector and (X^a - -^(«-i)a) T 
denotes its transpose, and we set X% = for i < 0. More precisely, given a 
time shift r G R, we are interested in the behaviour of Va,t(N,N t+ .). In 
essence, Va ; t(JV,]V t+ .) can be viewed as a multivariate cross-correlogram 
across scales A of N: it can be consistently measured from empirical data 
and its limiting behaviour as T — > oo plays a key tool in understanding the 
second-order structure of linear functions of N across scales A. 

Theorem 3. In the same setting as in Theorem 2, let (At)t>o an d ( t t)t>o 
be two families of real numbers such that At > 0. If At/T — > and 
tt/T —t-0 as T —7- oo, we have 

Va t ,t(X, X Tt+ .) — va t ,t t — > as T —?■ oo in L 2 (P) 

where 

VA, T = (l-§) + S+ / d S dt(l-^^) H >(s)S^) T + 

poc poo 

+ ds(l- i£+d)+V>( 5 )E + / ds (1 - ^) + S </>( S ) T , (8) 
or equivalently, 

v A , T = / (l-^^) + (ld5o(ds) + V(s)rfs)S(ld«5o(^) + ^(i) T ^) 1 

where XI is f/te diagonal matrix such that Yin = ((Id — K) _1 //J. anc ^ ^ e 
function ip is given by (7). 
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Remark 1. On can check that va )T = E((N& — Nq)(N&+ t — -^r) T ) where N 
is a counting process of the (unique in law) stationary multivariate Hawkes 
process on R associated to \x and cp. Thus, another way to obtain va, t is to 
compute E((Na — Nq)(Na +t — N T ) T ) in the stationary regime, as in [13, 12, 
14] by means of the Bartlett spectrum of N, see [4]. However, the stationary 
restriction is superfluous and moreover, only very specific parametric form 
of (fij like exponential functions enable to carry such computations. 

Remark 2. Obviously, we have va, t = vT . 

Remark 3. For fixed r G R, we retrieve the macroscopic variance of Theorem 
2 by letting A — > oo. More precisely, we have 

va,t ->rsr T as A ->■ oo 

where T = (Id — K) -1 , and the effect of r vanishes as A — > oo. For all 
r/Owe have 

va,t -> as A — > 0. 

This convergence simply expresses the fact that the two processes ./V and 
iV T+ . cannot jump at the same time, producing a flat autocorrelogram for 
sufficiently small sampling mesh A. 

In the same way as Corollary 1 is obtained from Theorem 2, we have the 
following refinement of Theorem 3. 

Corollary 2. Finally, in the same setting as in Theorem 3, assume more- 
over that (A2) holds. Define 

~ f N t - i(Id - K)~V if t>0, 
* \ if t < 0. 

Va t j(^,Xt t +-) - va TiTt — » as T — ^ oo in L 2 (P), 
where va, t 5«wen 6y (8). 

5 Application to financial statistics 

5.1 The macroscopic trace of microstructure noise 

Following [17], we introduce a univariate price process S = (St)t>o by setting 

S = JVi - JV 2 , 
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where (Ni, N2) is a Hawkes process in the case d = 2, with 



M,t = " + J( ,t) - s)dN 2)S , 
A 2ji = f + J (o t) - s)dN 1>s 



for some f 6 K4. and <p : R + — )• R + . With our notation, this corresponds to 
having fi = {y, v) and 

'0 ^ 
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If 93 = 0, we find back a compund Poisson process with intensity v and 
symmetric Bernoulli jumps. This corresponds to the simplest model for a 
random walk in continuous time, constrained to live on a lattice, the tick-grid 
in financial statistics accounting for the discreteness of price at fin scales. 
Microstructure noise corresponds to the property that an upward jump of 
S will be more likely followed by a downward jump and vice versa. This 
phenomenon lays its roots in microeconomic analysis of price manipulation 
of agents [10, 23, 24, 28, 29, 21]. In our simple phenomenological setting, 
it will be reproduced by the introduction of the kernel ip, as empirically 
demonstrated in [17]. The question we can now address is the macroscopic 
stability of the model. Do we retrieve a standard diffusion in the limit 
T — > 00 for an approriate scaling of S and how does the effect of tp influence 
the macroscopic volatility? By Theorems 1 and 2, we readily obtain an 
explicit anwser: 

Proposition 1 (Macroscopic trace of microstructure noise). Assume that 
IMIl 1 = fo° f{t)dt < 1. Then 

(T-V 2 S T v,v G [0, 1]) -»• (aW v ,v G [0, 1]) as T -> 00, 

in law for the Skorokod topology, where {Wv)v^[o,\\ ^ s a standard Brownian 
motion and 

er - 



(i-|M| L1 )(i + IM| L i) 2 - 

Remark 4. Note that if we take (p = 0, we retrieve the standard convergence 
of a compound Poisson process with symmetric jump to a Brownian motion. 

Remark 5. By assumption, < ||y||x,i < 1 and a closer inspection of the 
function 

x ~* a(y, xf = — — — - for x G [0, 1) 

(1 — x)(l + x) z 
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reveals an interesting feature: for small microstructure effect (namely if 
x = || 1,1 l ess than 1/3) the effect of microstructure tends to stabilise the 
macroscopic variance in the sense that a(u, x) 2 < a(u,0) 2 , whereas beyond 
a critical value x ~ 0.61, we have a(u, x) 2 > a(u, 0) 2 and even o~(v, x) 2 — > oo 
as x — > 1. 

5.2 Macroscopic correlations for bivariate assets 

We now turn to a bivariate price model S = (Si, S2) obtained from a Hawkes 
process in dimension d = 4, of the form 

(S 1 ,S 2 ) = (N 1 -N 2 ,N 3 -N 4 ) 

with fi = (//I, fi2, fJ-3, /m) such that 

/il = /.f 2 and fi 3 = /i 4 

together with 

/0 h 0\ 

/i 

¥ '~ g 

\o 9 0./ 

The upward jumps of Si excite the intensity of the upward jumps of S2 and, 
in a symmetric way, the downward jumps of <Si excite the downward jumps 
of S2 via the kernel g : M + — > R + . Likewise, the upward jumps of S2 excite 
the upward jumps of Si and the downward jumps of S2 excite the downward 
jumps of Si via the kernel h : R+ — > M+. For simplicity we ignore other 
cross terms that could produce microstructure noise within the inner jumps 
of Si and S2 ■ This relativeley simple dependence structure at a microscopic 
level enables to obtain a non-trivial form of the macroscopic correlation of 
the diffusion limits of Si and S2. 

Proposition 2. Assume that \\h\\£i ||ff||ii < 1- The 2- dimensional processes 
T" 1 / 2 {Si t T v , S2 t Tv) vG [ x j converge in law as T — > 00 for the Skorokod topol- 
ogy to 

(Xi\ = yg / V \' 2 Wi + vl ,2 \\h\\ L i W 2 \ 

\X 2 ) (1 - \\h\\ L 4g\\ L ,f/ 2 \v\' 2 \\g\\ Ll Wi + 4 /2 W2j 

with 

v\ = m + \\h\\ L ifj,3, v 2 = H3 + \\g\\ L ifii. (9) 
and where (Wi,W2) = (Wij, W / 2,t)iG[o,i] ^ s a standard Brownian motion. 



13 



The proof is a consequence of Theorems 1 and 2 and is given in appendix. 

Remark 6. The macroscopic correlation between S\ and 52 is thus equal to 
the cosine of the angle of the two vectors 

("J . "2 IHU 1 ) and (^llflU 1 )^ )• 

Obviously, it is always nonnegative and strictly less than 1 since the deter- 
minant v\^ 2 v^ 2 (l — || /i || £i 1 1 5 1 1 l 1 ) °f the two above vectors is positive unless 
the Hi are all 0. 

5.3 Lead-lag and Epps effect through the cross correlations 
across-scales 

We keep up with the model and the notation of Section 5.2 but we now 
study the quantities 

Va,t(S\, 5i )T+ .), Va,t(S2,S2, t +-), Va,t(Si,S2 )T +-)- 

In particular, the quantity Va,t(Si, S^t+O is a powerful tool for the statisti- 
cal study of lead-lag effects, i.e. the fact that the jumps of Si can anticipate 
on those of S2 and vice- versa, see for instance [1, 30]. The Epps effect - i.e. 
the stylised fact statement that the correlation between the increments of 
two assets vanishes at fine scales - can be tracked down likewise. Theorem 
3 enables to characterise in principle the limiting behaviour of these func- 
tionals. This is described in details in Proposition 3 below. 

We consider g and h as functions defined on R by setting g(t) = h(t) = 
is t < 0. Assume that \\h\\ L i \\g\\i,i < 1. Then the series 

n>l 

converges in L 1 (M, dt). If / is a function on R we define / by f(t) = f(—t). 
We have 

Proposition 3. Assume that \\h\\ii \\g\\Li < 1. Let (At)t>o and (tt)t>o be 
two families of real numbers such that At > 0. If Ay/T — > and tt/T — > 
as T —7- 00 we have 

Va t ,t(Si, Si >Tt+ .) — Ch(At, tt) — > as T — ^ 00 in L 2 (P), 

where 

Cn(A,r) = t n fe || \ M - 7A * (<5o + F + F + F * F) * (i^Sq + v 2 ,h* h) (r), 

ll 1 1 1 1 1 y 11 1 
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with 

7A(x) = (l-|x|/A) + . 

We also have 

Va t ,t(Si,S2,tt+-) - Ci2(A T ,T T ) — > as T —?■ oo in L 2 (P) 

and 

Va t ,t(S2, 5 , i jTt+ .) - Ci 2 (At, —tt) -> as T -»■ oo m L 2 (P), 

with 

C 12 (A, r) = ^i^ ii 2 I, I, 7A * (5 + F + F + F * F) * (i/^ + i^) (r). 

II II.L- 1 - NjII/j-L 

Remark 7 (The Epps effect). For / G L 1 we have 7a * / — ► pointwise as 
A —> 0. Therefore we obtain 

Ci 2 (A, r) -> as A ->■ for every r G R 

and this characterises the Epps effect. The same argument, together with 
7a(0) = 1, yields the convergence 

2 v\ 

Cn(A,r)^- l {r=0} asA^O. 

Remark 8 (The Lead-Lag effect). Following [17] and as a consensus in the 
literature, we say that a lead-lag effect is present between S\ and 5 2 if there 
exists A > and such that 

C 12 (A,r) /C 12 (A,-r). 

Therefore, an absence of any lead-lag effect is obtained if and only if the 
function (5o + F + F + F * Fj ★ (i/ 2 /t + z^i<?) is even. This is the case if 
v\g = z^ 2 /i. Now, let e > 0. If g = /i ★ 5 £ and /ii = /i3, then 

Ci 2 (A,-r) =Ci 2 (A,r + e) for every A>0,rGM. 

This particular choice for /i and g models in particular the property that S± 
acts on S 2 in the same manner as 5" 2 acts on Si with an extra temporal shift 
of e. Since we always have 

Urn Ci 2 (A,r) = 0, 

T— ¥±00 

there exists To such that Ci 2 (A, tq + e) ^ Ci 2 (A, To), or in other words, we 
have a lead-lag effect. 
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Remark 9 (Macroscopic correlations). Since 7a — > 1 as A — > oo, we obtain 
the convergence 

Cu(A,t) -> izpiCTi^ / (<5o + F + ^ + i ? *F)*M) + ^*/*) 

= i- || ft || fil| g | | £ i (l + 2 / ^ + (/ F) 2 ) + hf) = Var(Xi). 

Likewise, we have Ci2(A,r) — >■ Cov(Xi,X 2 ) as A — >• oo. 

Remark 10. Finally, note that if we use the convenient parametrisation 
h(t) = a\ exp(— fat) and g(t) = a 2 exp(— fat), assuming further \\h\\ L i — 
0:102/ (/Si A) < 1) then standard computations yield the explicit form 

F(t) = (exp(-p 2 t) - exp(-pii))l]R + (t) 

Pi ~ P2 

with 

Pi = \ (pi + P2 + ~ P2Y + 4ai« 2 ) , 

P2 = r (ft + & - V(A - ft) 2 + 4«ia 2 ) 

and we have < p 2 < Pi- This allows to obtain a close formula for Cu and 
Ci 2 . We do not pursue these computations here. 

6 Preparation for the proofs 

In the sequel, we work in the setting of Sections 2 and 3 under Assumption 
(Al). 

Lemma 2. For all finite stopping time S one has: 

E(N S ) = pE{S) + e(J <p(S - t)N t dt) (10) 

E(N S ) < (Id - K)~ 1 /j,E(S) componentwise. (11) 

Proof. Recall that (T p ) p >i denote the successive jump times of N and set 
S p = S AT p . Since the stochastic intensities Aj are given by (4) one has 

E(N Sp ) = pE(S p ) + E( [ SP dt [ cp(t - s)dN s ) . 

V0 J(0,t) J 
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Moreover, by Pubini theorem 

rS p p p pS p 

/ dt (p(t-a)dN s = / (/ ip(t - s)dt)dN s 
Jo J(o,t) J[o,S p ) Js 

p pS p -s 

= (/ <p{t)dt)dN a . 

J[0,S p ) Jo 

Now, integrating by part with <&(i) = J * (p(s)ds, we can write 

= &(0)N Sp - $(S P )N = [ <f>(S p - t)dN t - <p(S p - t)N t dt. 

J(o,s p ] Jo 

Remark that both sides of the above equality are finite since Yli=i Ni,s p < P- 
We obtain 

E(N Sp ) = fiE(Sp) + E(JJ f(S p - t)N t dt) 
and derive (10) using that N$ t N$ as p — > oo and 

/ * <p(S p -t)N t dt = f ' ' <p(t)N Sp - t dt t / tp(t)N S -tdt = f <p(S-t)N t dt. 
Jo Jo Jo Jo 

We next prove (11). We have 

E(N Sp ) = e(s p h + J J cp(S p - t)N t dt) 

< E(S p )fi + e(^J f(S p — t)N t dt) componentwise, 
= E(S p )fi + KE(N Sp ). 

By induction 

E(N Sp ) < (ld + K + -.-+K n - 1 )E(S p )» + K n E(N Sp ) 

componentwise for all integer n. On the one hand, Ylt=i ^i,S p < P- On the 
other hand, since p(K) < 1 we have K n — > as n — > oo and 

oo 

^K" = (Id-K)- 1 , 

n=0 

therefore 

E(N Sp ) < (Id-K)- 1 ^)^ 
This readily yields (11) since E(N S ) = \im p E(N Sp ). □ 
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Let ip n , n > 1, and ijj = ^ n>1 i/J„ be defined as in Theorem 2. By- 
induction it is easily shown that J °° (p n (t)dt = K n for all n. Therefore 
J °° tj}(t)dt = ^2 n>1 K n is finite componentwise by assumption (Al). We 
next state a multivariate version of the well known renewal equation, which 
proof we recall for sake of completeness. 

Lemma 3. Let h be a Borel and locally bounded function from R+ to M. d . 
Then there exists a unique locally bounded function f : — > M, d solution to 



given by 



f(t) = h{t) + I cp(t- s)f(s)ds Vi > 0, (12) 
o 



f h {t) = h(t) + rjj{t- s)h(s)ds. 



o 



Proof. Since tp E L l (dt) and h is locally bounded, the function fh is locally 
bounded. Moreover fh satisfies (12). It follows that 

cp(t — s)fh(s)ds = / (f(t — s)h(s)ds + / dsip(t — s) / ip(s — r)h{r)dr 



o 



/ <p(t — s)h(s)ds + dr dscp(t — s)if>(s — r)h(r) 
Jo Jo Jr 



"0(t — r)h{r)dr 



o 



since f Q ip(t — s)tp(s)ds = tp(t) — y>(t). As for the uniqueness, if / satisfies 
(12) then 



fh(t) ~ f(t)= / <p(t-s)(f h (s)-f(s))ds 
Jo 

thus if gi(t) = \f h ,i(t) - fi(t)\, 1 < i < d, one has 

g(t) < / f{t — s)g(s)ds componentwise, 
Jo 

which yields 

poo poo 

/ g{t)dt < K / g(t)dt componentwise. 
Jo Jo 

Since p(K) < 1 it follows that f = fh almost everywhere. Therefore 

/ (f(t — s)f(s)ds = / (p(t — s)fh(s)ds for all t 
Jo Jo 

and thus f = fh since both function satisfies (12). □ 



18 



Define the d-dimensional martingale (Mt)t>o by 

M t = N t - f X s ds with A = (Ai,...,A d ). 
J o 

Lemma 4. For all t > 0: 

E{N t )=t/i+(f ip(t-s)sds)n, (13) 
Jo 

N t - E(N t ) = M t + [ ij>(t - s)M s ds. (14) 
Jo 

Proof. By (10) of Lemma 2 and Fubini theorem, we get 

E(N t ) = tfi+ [ tp(t - s)E(N s ) ds. 
Jo 

Besides, t E(Nt) is locally bounded in view of (11). Applying Lemma 
3 we obtain (13). The second formula follows from Lemma 3 and the fact 
that, if X t = N t — E(N t ), representation (13) entails 

X t = M t + [ ip(t - s)X s ds. 
Jo 

□ 

7 Proof of Theorem 1 

Lemma 5. Let p € [0,1) and assume that J °° t p (p(t) dt < oo component- 
wise. Then 

1. If p < 1, we have 

Tp(t- x E{N Tv ) - v{ld - K)~V) ^0 as T -> oo 

uniformly in v £ [0, 1]. 

2. If p = 1, we have 

r(i#(iv T )-(id-K)-V) 

-> - (Id - K)" 1 ( / ttp(t) dt) (Id - K)~V as T — ^ oo. 
Jo 
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We first prove that L°° t p ij)(t) dt < oo componentwise. For n > 1, setting 



Proof. Let p G [0, 1] and assume that t p ip{t) dt < oo componentwise. 
We first prove that J °° t p if)(t) dt < 
A n = Jq 00 t p (f n (t) dt, we can write 

A n+ i = J cp(t-s)ip n (s)d s yt 



oo , roo 



{t + s) p ip(t)dt)v n (s)ds 



o vo 

oo 



/*oo poc 

< / t p (p(t) dt K n + K / s*V n (s) (is with equality if p = 1, 
Jo Jo 



= AiK n + KA n . 

Therefore for all integer N, 



N N-l N-l 



and 



Y A„ < Ai + Ai Y K n + K Y An, 

n=l n=l n=l 

AT-1 7V-1 

(Id - K) Y A n + A N <A l + A x Y K r 

n=l n=l 
AT-1 JV-1 



^ A n < (Id - K)- 1 (A 1 + Ai Y K " 

n=l n=l 

Letting iV-^oowe derive 



t^^(t) dt = ^ A n < (Id - K)- 1 Ai(Id - K)" 1 

n>l 

with equality if p = 1. From (13) it follows that for all t> G [0, 1]: 

u(Id - K)"V - ^E{N Tv ) = [v ^(t)dt + - J t*/>{t)dty (15) 
Since t v il){t) is integrable, we have 

/"OO poo 

T p / ^(t)dt < v l ~ p / f^(t) dt as T -> oo 

and this convergence is uniform in u G [0, 1] in the case p < 1. It remains to 
prove that if p < 1 , we have 

cT 



T l ~ p J 



[ ttp(t)dt ^ as T -)• oo. 
Jo 
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With G(i) = JjJ s p ip(s)ds, integrating by part, we obtain 

T 1-p G(T) = f tip{t) dt + (l-p) [ t~ p G(t) dt 
Jo Jo 

and 

-L £ ^(t)tdt = G(T)-^£ r*G(t) dt. 

Since G(t) is convergent as t — > oo we finally derive that the right-hand-side 
in the above equality converges to as T — > oo. □ 

Denote by || • || the Euclidean norm either on R d or on the set of d x d 
matrices. 

Lemma 6. There exists a constant such that for all t, A > 0: 
E( sup \\M S -M t f) <C^A. 

t<s<t+A 

Proof. Doob's inequality yields 

d 

E( sup \\M S -M t f) <4V J B((M M+A -M M ) 2 ). 

t<s<t+A ~^ 

For each i G {1, . . . , d}, the quadratic variation of the martingale (Mij)t>o 
is 

[M u M t ] t = £(M M - M M _) 2 = N i>t . 

s<t 

Thus we have 

E((Mi, t+ A - Mi, t f) = E(Ni, t +A - Ni, t ). 

Besides, in view of Lemma 4 and the fact that J °° tp(t)dt = (Id — K) _1 — Id, 
we obtain 

E{N t+A - N t ) < A(Id - K)~V componentwise. 

□ 

Completion of proof of Theorem 1. Lemma 5 with p = implies that it is 
enough to prove the following convergence 

T _1 sup \\N Tv - E(N Tv )\\ ->0 as T ^ oo (16) 
ve[o,i] 
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almost surely and in L 2 (P). Thanks to (14) of Lemma 4, we have 
11^-^(^)11 < (1+ / 

^e[o,i] 



sup \\N Tv - E(N Tv )\\ < (1 + / ||^(t)||dt)sup||M t ||, 
wero.il Jo t<T 



< C v sup 1 1 M<| 

t<T 



since i/> is integrable. Moreover 



£(sup||Mi|| 2 ) < C^T 

t<T 

by Lemma 6, therefore convergence (16) holds in L 2 (P). In order to prove 
the almost-sure convergence, it is enough to show that 

T" 1 sup H-Mr^H —t-0 as T — > oo almost-surely. 

«6[0,1] 

Let 

Z t = (Z 1>t ,...,Z d>t ) = [ -^—dM s . 
The quadratic variation of the martingale Z^ satisfies 

[z i} Zi] t = ( z i,s - z t , s -f = ^ (7^i)2 dNi 



0<s<t 



and moreover, using integration by part and (13), we have 

f 1 f N 

V(o,oc) (a + 1) ; V ( o,oo) (1 + s) 3 ; 

Therefore limt^oo Zt exists and is finite almost surely. It follows that 



1 1 



t+l t+1 



Mt = Zt / Z s ds — > as T — > oo almost surely. 



o 



We deduce that almost surely, for all family vt G [0,1], T > 0, such 
that Tvt — > oo the convergence Mt Vt /T — > holds. Moreover, we have 
Mtv t /T — > if swp T TvT < oo. In other words T~ 1 Mt v uniformly 
uniformly m v G [0,1], almost-surely. The proof of Theorem 1 is com- 
plete. □ 
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8 Proof of Theorem 2 

Let W = (Wi, . . . , Wd) be a standard d-dimensional Brownian motion. For 
i = 1, . . . , d, put ai = (Sjj) 1/2 - 

Lemma 7. The martingales := (T -1 ' 2 Mtv)v€\q,i] converge in law for 

the Skorokod topology to (ciWi, . . . , cr^W^). 

Proof. According to Theorem VIII-3.11 of [19], since the martingales 
have uniformly bounded jumps, a necessary and sufficient condition to ob- 
tain the lemma is: for all v £ [0, 1], for all 1 < i < j < d 

[m\ t \mI t \ -> afv, [mP,mP] v ^0, as T ->■ oo in probability. 
We have 

[MP,MP} V = ^N i>Tv -> afv in L 2 (P) by Theorem 1 

and 

[MP,Mj% = 

since the processes iVj for 1 < i < d, have no common jump by construction. 

□ 

Completion of proof of Theorem 2. Set 

XP=T-V*{N TV -E{N TV )). 
In view of Lemma 7, it is enough to prove that 

sup \\X^ - (Id - K) _1 M„ (T) || ->• as T -)• oo in probability. 

ve[o,i] 

By Lemma 4, we have 

Zf ) = Mf ) + f T^(Tu)mZ{ du, 
Jo 

hence we need to prove that 

sup || / T^{Tu)M^J u du - ( / ^(t)^)M^ T )|| -> as T -> oo (17) 

t>€[0,l] •/ J 
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in probability. We plan to use the fact that ip is integrable and the C- 
tightness of the family (M^)t>o- The tightness is a consequence of Lemma 
7 and reads: 

Ve>0 limsupP( sup || M^ T) - M*p \\ > e) -> as 77 -> 0. (18) 

(T) 

using also that Mq =0. For 77 > and u 6 [0, 1] we have 

I' 1 

T^(Tu)M { Jj u du\\ < sup HHP]] / r|^(Tu)||dti 

tiA'jy 0<u<l Jr\ 

poo 

< sup \\MPW / ||^(t)||dt^0 

0<m<1 JTr/ 

as T — >• 00 in probability, since sup 0<u<1 \\Mu \\ is bounded in probability 
and ->■ as T -»■ 00. Moreover 

pvAr/ roo 

11/ ^(^(M^-M^J^H < sup \\MP-M { J ] \\ / ||V(t)||dt, 

io \u-u'\<r) JO 

therefore, in order to prove (17) it suffices to show that for all 77 > 0, the 
convergence 

/•oo i-vArj 

sup ||(/ ip(t)dt- Til>(Tu)du)Mjp\\ — > as T — > 00 
ug[o,i] Jo Jo 

holds in probability It readily follows from (18) and the upper-bound 

rvArj 



a OO pV/\T] . 

ij)(t)dt - J T^(Tu)duj Mjp I 

\mt)\\dt sup||Mf)|| ifu>a 



< < 

/>oo 

/ \mt)\\dt sup||Mf)|| if 

, JO u<<5 

with < (5 < r). □ 

Proof of Corollary 1. By (A2), Lemma 5 with p = 1/2 yields 
r V2( r -i£(7v- Tv ) _ v ( Id _ K)~V) as T — > 00 

uniformly in u G [0, 1]. Moreover, by (Al), Theorem 2 yields 

T 1/2 (T^N Tv - T~ l E{N Tv )) -»• (Id - K)" 1 ^ 2 !^ 

in distribution as T — >• 00 and the result follows. □ 
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9 Proof of Theorem 3 

Set 

Y t = ! ^{t-s)M s ds (19) 
J o 

in order that X = M + Y , see Lemma 4. For all < e < n < 1, for all 
integer 1 < ko < T/A, define 

^k ,e,n( X )T,A = ^ ^ (^(fe-l)A+r/A - ^(As-1)A+ea) ■ 
k=ko 

Lemma 8. There exists a function T £n tp(T) such that £^ i¥ ,(T) — ^ as 
T — > oo and swc/i i/iai /or all < e < rj < 1 anci all integer 1 < ko < T/ A, 
we have 

E{\\(1 - ^)(rj - e)(Id - K)"V - ZW^)r,A|| 2 ) < WT). 
Proof. First we prove that 

E(||P fc0ie>r? (X) T , A || 2 ) < C^T- 1 (20) 
for some constant that depends on fj, and <£> only. Using 

^(fc-l)A+??A - ^(jfc-l)A+eA = M (k~l)A+rjA ~ -^(fc-l)A+eA 
roc 

+ / dsip(s)(M( k _ 1)A+vA _ s - M {k _ 1)A+eA _ s ) 
Jo 

and the fact that ip is integrable, it suffices to prove that 
£(||l? feo>E ^(M._ s ) T)A || 2 ) < C^T-\ 

Set 

n 

Sn = ^(^(fc-l)A+r;A-s ~ -^(fc-l) A+eA-s) • 
k=l 

Clearly (S n ) n >i is a discrete martingale thus 

LTAJ 

^(II^AJ -5fc _i|| ) = E(\\M( k _i) A+vA _ s - M^-^A+eA-all )■ 
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By Lemma 6 we have E (US' [TAJ || 2 ) < C^,ip(j} — £)T an d (20) follows. It 
remains to prove that 

||(1 - ^)(r? - e)(Id - K)~V - E(V k0>E>ri (N) TA ) \\ < g^T) 

where £u,ip(T) — > as T — > oo. By Lemma 4, we have 

E{N(k-l)A+r)A - N (k-l)A+eA) 
r(k-l)A+eA 

= (t] — e)A(h+ / dn/>(r)//J 



+ / drip(r)((k - 1)A + r]A - r)fi 

J(k-l)A+eA 

POO 

(rj - e)A(Id - K) _1 /x - (v ~ / (irV'(r) 

J(fc-l)A+eA 

p(k—X)A+7]A 

+ / drip(r)((k-l)A + r]A-r)n. 

J(k-l)A+sA 



Finally 

r(k— l)A+r;A >>oo 

/ dri/>(r)((k — 1)A + ??A — r) < (77 — e)A / drip(r). 

J(k-l)A+eA J(k-l)A+eA 

We conclude noting that 

A L^/aj A |T/AJ 

— V / dri/>(r) = — / dr^(r) V l { ( fc _i)A+ e A<r-} 

J fc=fco ./(fc-l)A+eA ./(fco-lJA+eA k=ko 

<-/ dr-0(r)(rAT), 



T . 

This last quantity converges to as T — > 00 by an argument similar to the 
end of proof of Lemma, 5, using that ip is integrable. □ 

Lemma 9. There exists a constant C^ tV such that for all t, h > 0, we have 
E( sup \\M S - M t \\ 4 ) < (h + h 2 ) 

t<s<t+h 

Proof. According to the Burkholder-Davis-Gundy inequality, we have 

d 

E( sup \\M S - M t \\ A I & t ) < C Y^E({[M u Mi] t + h - [M^M^f \ & t ) 

t<s<t+h ■ ; 

— — 2 = 1 

d 



i=l 
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hence 

E( sup \\M S -M t t) <C{E{\\X t+h -X t \\ 2 ) + \\E{N t+h -N t )\\ 2 ). 

t<s<t+h 

By Lemma 4, we have 

E^Nt+h — Nt) < h(Id — K)" 1 ^ componentwise 

and 

(E(\\Xt+ h - Xtf)) 1 ' 2 < U{s)\\ (E(\\M t+h _ s - M t || 2 )) V2 ds. 

The conclusion follows using Lemma 6 and the fact that ift is integrable. □ 

With the notation introduced in Section 3, the quantity Va,t(-M- 3 , M._ t+ T ) 
is equal to 

j \T/A\ 

f ^2 { M kA-s ~ M (k-l)A-s) {MkA-t+r ~ M(k-l)A-t+r) ■ 
k=l 

Lemma 10. For all s, t > we have 

V AT , T (M._ s ,M._ t+TT ) - (1 - ^^) + E -> as T — ^ oo m I 2 (P) 
and 

^(||V ATiT (M_ s ,M_ i+TT )|| 2 ) < C^. 

Proof. We start with proving preliminary estimates. Let &i, &2> ^4) ^2 be real 
numbers such that 61 < 62, 6' x < 6 2 an d (£>i, fo) D (6^, 6 2 ) = 0- Using that M 
is a martingale, we successively obtain 

LT/Aj 

S (ll r 1 X] ( M (fc-l)A+fe 2 ~ M (fc-l)A+fei)( M (fc-l)A+fe' 2 ~ M (fc-1)A+6'J || ) 
fe=l 

j LT/Aj 

= ^2 ^(11 ( M (fc-l)A+fe 2 ~ ^(fe-l)A+6i) {M(k-l)A+b' 2 ~ ^(fc-l)A+6i) || ) 

fc=l 
j L^/AJ 

= f2 - E '(ll M (fc-l)A+62 ~ -^(fc-l)A+6i|| \\{ M (k-l)A+b' 2 - M (k-l)A+b'l) || )• 

fc=l 
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By Cauchy-Schwarz, this last quantity is less than 
1 LT/AJ 

if2 E ( E (U^i)A+b 2 -M {k _ 1)A+b J 4 )E(\\(M { ^^ 

k=l 

which in turn, using Lemma 9, is bounded by 

C^^(6 2 - h + (b 2 - hf) 1 ' 2 ^ - b[ + (b' 2 - &;) 2 ) 1/2 . (21) 

Moreover, if b 2 — b\ < A, using that [M, M] = diag(iV) and Lemma 9, we 
obtain that 



E 



1 l T / A \ 

— dia g( N (k-l)A+b 2 - iV(fc-l)A T +6i) 

fc=l 



LT/Aj 



~T ^ ( M ( fc " 1 ) A + b 2 _ Af (fe-l)A+6i)( A ^(fe-l)A+6 2 ~ M (k-l)A+b 1 ) T ) 
fc=l 

is less than 

C^^-(6 2 -6i + (6 2 -6i) 2 ). (22) 

We are ready to prove Lemma 10. It is a consequence of Lemma 8 and the 
fact that there exists a 2 = a 2 (s, t, At, tt) < 03 = 0,3(3, t, At, tt) such that 

03 - «2 = {At — \t — s — t t \) + , as T ~^ 00 

holds, together with the estimate 

j Lt/a t j 
e (\\t S dia s(% 

fc=l 

<C^i(l + A T ). 



(23) 



Indeed, if t — tt + A^ < s or s < t - rj - At then the upper bound we 
obtained in (21) entails 

£(||v At , t (m_ s ,m._ 4+7T )|| 2 ) <c^1(i + a t ). 

Let us first consider the case £ — tt < s < t — tt + At- Set 

ai := s < a 2 := — i + tt < £13 := A^ — s < 04 := A^ — t + ry. 
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We use the following decomposition 

(M fc A T - s - M (k _ l)AT _ s ) (M fc A T -t+T T — M( fe _l)A T -t+T T ) 
= (%-l)A T +a 3 ~~ -^(fc-l)A T +a 2 ) { M {k-l)A T +a 3 ~ M {k-l)A T +a 2 ) T 
+ {M(k-l)A T +a 3 - M{k-l)A T + ai ) { M (k-l)A T +a 4 ~ M (k-l)A T +a 3 ) T 
+ {M(k-l)A T +a 2 - -^(fc-l)A T +oi) (-^(fc-l)A T +03 _ -^(fc-l)A T +a 2 ) T - 

On the one hand, (21) readily yields that both 
1 Lr/A T j 

E (\\f ^2 ( M (fc-l)A T +a 3 - M (fc-l)A T +ai)( M (fc-l)A T +a4- M (fc-l)A T +a 3 ) || ) 
fc=l 

and 

1 Lr/A T j 

E (\\f ^ ( M (k-l)A T +a 2 - M (k-l)A T +a 1 ){M(k-l)A T +a3- M (k-l)A T +a 2 ) || ) 
k=l 

are less than (7^^,(1+ At) /T. On the other hand, by (22), the same estimate 
holds for 

1 IT/A T \ 

E (\\t Yl dia g(%-l)A T +a 3 - %-l)A T +a 2 ) 
fc=l 

1 LT/ATJ T 2 

~^ X] ( M (fc-l)A T +a3- M (fc-l)A T +a 2 )( M (fc-l)A T +a 3 - M (fc-l)A T +a 2 ) || ), 
k=l 

therefore (23) holds in that case. If now t — tt — At < s < t — tt, setting 

di := — t + tt < a2 := — s < «3 := Ay — t + tt < 04 := A^ — s, 
one readily checks that (23) holds using the same arguments. □ 

Completion of proof of Theorem 3. Since X = M + Y by (19), we have the 
following decomposition 

V AjT (X,X T+ .) 

= V A ,T(M,M T+ .)+V AiT (Y,Y T+ .) + V At T(M,Y T+ .)+V A , T (Y,M T+ .). 
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Setting M t = for t < 0, we can write Y t = J* °° -0(s)Mt_ s ds, therefore 

(YiA - ^(i-i)a)(^A+t - Y(i-l)A+r) T 

= f dsdt^(s)(M iA _ s - M {l _ 1)A _ s ) (MiA-t+r - M {i _ 1)A _ t+T ) T ^(t) T , 
Jr 2 



hence 



Likewise 



V a ,t(Y,Y t+ .) = [ dsdt^(s)V AjT (M._ s ,M.^ +T )iP(t) T . 

JR 2 + 

/•oo 

V A ,T(M,y r+ .)= / dtV AtT (M,M.^ t+T )iP(t) T , 
Jo 



oo 



V AiT (y,M T+ .)= / dt*l>(t)V A , T (M._ t ,M. +T ) 



o 



In view of Lemma 10 and the fact that ip is integrable, by Lebesgue domi- 
nated convergence theorem, we successively obtain 

V At , t (Y,Y Tt+ .)- [ dsdt^(s)(l - \t - s - t t \/A t ) + S^(t) T -> 0, 

Jr 2 + 

V AtiT (M,7 Tt+ .)- / ^(l-|t-r T |/A T ) + 5]^(t) T ^0, 
Jo 

and 

/•oo 

V At ,t(^M tt+ .)- / dt^(t)(l-\t + T T \/A T ) + V^0 
Jo 

in L 2 (P) as T — 7- oo. The proof is complete. □ 
Proof of Corollary 2. In view of Theorem 3, we need to show 

V At ,t(X,X tt+ .) - V At , t (X,X tt+ .) -> as T oo in L 2 (P). 
By Cauchy-Schwarz inequality, this convergence is a consequence of 

sup- V (r ?i (r + A;A T )-r ?i (T + (A;-l)A T )) -> 0, 1 < i < d, 

where 77(f) = (t(Id - K)~V - l{t>o}- By (15), for t > 0, we have 

the decomposition 

= (t J tp(s)ds + J ip(s)sdsjn 
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therefore the function r\ is absolutely continuous and we have 

/oo 
if>(s)ds[il {t > Q} . 

We derive 

y E (W + fcA T )-^(r + (fc-l)A T )) ^( S ) 2 d,. 
fc=i ^° 

It remains is to prove 

poo 

/ ||V(t)|| 2 <«<oo. 

■/ 

We have 

/oo 
s^ 2 iP(s)dsfi and h'WII <^~ 1/2 

since s ^ s 1 ^ 2 ip(s) is integrable. Finally 

h'(t)\\ 2 dt<c^ dtr 1 ' 2 dsU( s )\\ 

Jo Jt 

poo 

= 2C^ cZs||V(s)||s 1/2 < oo 

J 

and the result follows. □ 

Appendix 

Proof of Proposition 2 

/ /ft \ 

The spectral radius of K = I ? q q ^ J is equal to / g J h, and we have 

V o f g o o / 

/ l o fh o \ 



therefore 



( H - K >"'-Trk(i) 



31 



where v\ and v 2 are given by (9). Set X = N — E{N). By symmetry, 
we have N = {N 2 ,N 1 ,N 4: ,N 3 ) in distribution, thus E(N 1>t ) = E(N 2)t ) and 
E(N 3)t ) = E(N 4)t ) for all t. Consequently 

Si = X 1 - X 2 and S 2 = X 3 - X 4 . (24) 

According to Theorem 2 the processes T~ 1 / 2 Xt v converge in law to the 
process Y v = (Id - K.y l Y. l l 2 W v with 

-I / v\ \ 

v x I 

i/ 2 o I • 

v 2 I 



l-fhfg 



Therefore the processes T 1 / 2 (Si > t v , S 2) tv) v ^ ^ converge in distribution to 
m — ^2)^3 — YA and Proposition 2 is proved. 



Proof of Proposition 3 

From (24) it follows that 

V At ,t(Si,S 1iTt+ .) = (i-ioo)V At ,t(X,X, 
and 

Va t ,t{Sx,S 2 , Tt+ .) = (i -ioo)V At ,t(X,X tt+ .) 
Consequently Proposition 3 follows from Theorem 3 with 



Cn(A,r) = (i -i o o) v A)5 



o / ' 





C 13 (A,r) 



i-ioo)v Aj5 I 



It remains to compute C\\ and C\ 2 . First we compute tp = Yln>i ^Pn- We 
readily check that for all s,t > 0: ip{t — s)(p(s) = h(t — s)g(s)Id. Thus 
i~p 2 = (h-kg)Id. We derive 

(h±g)* n *h o 

(h*g)* n ±h 

(h*g)* n *g 

{h*g)* n *g 



t P2n = i h *9Y ni L&i <P 2 n+l 



and we obtain 



n>l n>0 

/ F (5 + F)*h 

F (5 +F)*h 

(5 + F)*g F 

V (5o + F)*5 F 



\ 
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where F = Y2n>i(h * d)*" '■ Set F(ds) = So(ds) + F(t)ds. Standard compu- 
tations yield 

/ F(ds) ^0 F*h(s)ds \ 

F(ds) ^0 F*h(s)ds 

F*g(s)ds F(ds) ^0 

\ F*g(s)ds F(ds) J 

and 



ld5 (ds) + ij)(s)ds 



(idJo(ds) + V(s)ds)E(ld«J (dt) +tp{t)*dt) = 

/an(ds,dt) a^ds, e?t) \ 

1 an(ds,dt) ai 3 (ds,dt) 

1 — fhfg a 3 i(ds,dt) a 33 (ds,dt) 

\ a 3 i(ds,dt) a 33 (ds,dt)J 



with: 



Therefore 



and 



a u (ds,dt) = v l F(ds)F{dt) + v 2 F * h{ds)F * h{dt), 
a 13 (ds, dt) = * #(dt) + ia^ * h(ds)F(dt), 

a 31 (ds, dt) = v 2 F{ds)F * /i(cft) + z^F * g(ds)F(dt), 
a 33 (ds, dt) = v 2 F(ds)F(dt) + z^F * * g(dt). 

Cii(A,r)=/ 7 A (i-s-r)2aii(ds,^), 

J[0,oo) 2 

C 12 (A,r)=/ 7 A (t-s-r)2a 3 i(cZs,cZt). 

J[0,oo) 2 



To complete the proof of Proposition 3, it suffices to use that for two finite 
measures \i and v on R one has (for all A > 0, r G R) 



7a(* — s — r)[i(ds)v(dt) = 7a * f * /2( r ) 
where /I is the image of /i by x ~^ — x. 
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